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We calculate glueball and torelon masses as well as the 
lowest lying hybrid potential in addition to the static ground 
state potential in lattice simulations of QCD with two flavours 
of dynamical Wilson fermions. The results are obtained on 
lattices with 16 3 x 32 and 24 3 x 40 sites at j3 = 5.6, corre- 
sponding to a lattice spacing, a -1 = 2.65^g GeV, as deter- 
mined from the Sommer force radius, at physical sea quark 
mass. The range spanned in the present study of five different 
quark masses is reflected in the ratios, 0.83 > m-n /m p > 0.57. 

PACS numbers: ll.15.Ha, 12.38.Gc, 12.39.Mk, 12.39.Pn 



I. INTRODUCTION 

The expectation that gluons form bound states, so- 
called glueballs, is as old as QCD itself jjj. Indeed, such 
states make up the spectrum of pure SU(3) gauge theory 
and have rather precisely been determined in quenched 
lattice simulations ^-|J. However, an unambiguous ex- 
perimental confirmation of their existence is still miss- 
ing. Lattice simulations have provided an explanation for 
this situation: most of the glueballs come out to be quite 
heavy (2|4j and will therefore manifest themselves as very 
broad resonances with many decay channels. This holds 
in particular for the most interesting glueballs, those with 
spin-exotic, i.e. quark model forbidden, quantum num- 
bers. 

However, a conservative interpretation of results ob- 
tained in the valence quark (or quenched) approxima- 
tion to QCD urges us to expect (see e.g. Rcf. ||) a 
predominantly gluonic bound state with scalar quan- 
tum numbers, J PC = ++ , and mass between 1.4 GeV 
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and 1.8 GeV. Indeed, in this region more scalar reso- 
nances have been established experimentally than a stan- 
dard quark model classification would suggest [pj . When 
switching on light quark flavours the difference between 
a glueball that contains sea quarks and a flavour singlet 
(isoscalar) meson that contains "glue" , sharing the same 
quantum numbers, becomes ill defined. In general such 
hypothetical, pure states will mix with each other to yield 
the observed hadron spectrum. This situation is different 
from that of flavour non-singlet hadrons which, if we ig- 
nore weak interactions, have a well defined valence quark 
content. Such "standard" hadrons might become unsta- 
ble but otherwise retain most of their qualitative proper- 
ties when sea quarks are included. Indeed ratios between 
masses of light standard hadrons, calculated within the 
quenched approximation, have been found to differ by 
less than about [Q 10 % from (full QCD) experiment. 

One can utilise quenched lattice data on the scalar 
glueball and isoscalar ss as well as uu + dd mesons as 
an input into phenomenologically motivated mixing mod- 
els [|S|. However, such models can only provide a crude 
scenario for the actual situation that will be clarified once 
the entire spectrum of full QCD has been calculated. 
A first attempt in this direction has been made by the 
HEMCGC Collaboration a decade ago |9j and recently 
first results from the UKQCD Collaboration on flavour 
singlet scalar meson/glueball masses in QCD with two 
light sea quark flavours have been reported |l(| . 

At present the main target of simulations involving 
sea quarks is to pin down differences with respect to the 
quenched approximation and to develop the methodol- 
ogy required when a decrease of sea quark masses below 
m^/rrip ~ 0.5 will become possible with the advent of 
the next generation of supercomputers. One would ex- 
pect sea quarks not only to be important for glueball- 
meson mixing and for an understanding of the spectrum 
of isoscalars but for flavour singlet phenomenology in gen- 
eralOO. 
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The potential between static colour sources at a sepa- 
ration, r, is one of the most precisely determined quan- 
tities in quenched lattice studies ]l^-[l5|. Two physi- 
cal phenomena are expected to occur when including sea 
quarks, one at large and the other at small distances. The 
former effect is known as "string breaking" (see Ref. flif l 
and references therein): once r exceeds a critical value, 
creation of a light quark- antiquark pair out of the vacuum 
becomes energetically preferable and the QCD "string" 
will break; the static-static state will decay into a pair of 
static-light mesonic bound states, resulting in the com- 
plete screening of colour charges at large distances that 
is observed in nature. The potential will approach a 
constant value at infinite r. On the other hand, the 
presence of sea quarks slows down the "running" of the 
QCD coupling as a function of the scale with respect to 
the quenched approximation: when running the coupling 
from an infra-red hadronic reference scale down to short 
distances, the effective Coulomb coupling in presence of 
sea quarks should, therefore, remain stronger than in the 
quenched case. 

Exploratory studies |l7|]l§| ] of the zero temperature 
static potential with dynamical fermions focussed on the 
question of string breaking. However, no statistically sig- 
nificant signal of colour screening has been detected. A 
difference between the potentials with and without sea 
quarks at short distances has then been reported by the 
SESAM Collaboration in Refs. Later studies by 

the UKQCD, CP-PACS and M1LC collaborations qual- 
itatively confirm these findings |^2|-p4|, employing dif- 
ferent lattice actions and, in the latter case, a different 
number of flavours. 

In this paper we present a combined analysis of 
SESAM and T^L @ data (as obtained on 16 3 x 32 and 
24 3 x 40 lattices, respectively) on the static potential, 
torelon and glueball masses. Moreover, for the first time, 
we have determined the II U hybrid potential [ p6| in a 
simulation including dynamical sea quarks. Preliminary 
versions of the results presented here have been published 
in Refs. @,|^@^^]Jl|2|§. 

The article is organised as follows: in Sect. || we sum- 
marise simulation details and our parameter values and 
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describe the measurement techniques applied. In Sect 
we discuss how the glueball and torelon masses and the 
potentials are determined. We also present the quality 
of the ground state overlaps achieved and the forms of 
the respective creation operators that were found to be 
most suitable. Subsequently, we present our physical re- 
sults on the potentials, torelons and glueballs in Sect. [j"v| , 
before we conclude. 



II. LATTICE TECHNOLOGY 

We analyse samples of configurations that have been 
generated by means of the hybrid Monte Carlo (HMC) 
algorithm using the Wilson fermionic and gluonic actions 



with nj — 2 mass degenerate quark flavours at the inverse 
lattice coupling, (3 = 5.6. This is done on L\L T = 16 3 32 
as well as on 24 3 40 lattices at the mass parameter values, 
K = 0.156,0.1565,0.157,0.1575 and 0.158. The corre- 
sponding chiralities can be quantified in terms of the ra- 
tios HJ, m w /TO p = 0.834(3), 0.813(9), 0.763(6), 0.704(5) 
and 0.574(13), respectively. The simulation parameters 
are displayed in Tab. |, where the Sommer scale [j3l| , 
tq w 0.5 fm, is defined through the static potential, V(r), 



dV 
dr 



1.65. 



(1) 



r=r 



At each k value 4,000-5,000 thermalised HMC tra- 
jectories have been generated. Great care has been 
spent |52 33 on investigating autocorrelations between 
successive HMC trajectories for the topological charge 
as well as for other quantities including smeared Wil- 
son loops and Polyakov lines, similar to those used in 
the present study. The operators that enter the glueball 
and torelon analysis have been separated by two (eight at 
K = 0.1565, five at k — 0.157) HMC trajectories; smeared 
Wilson loops that are required for the potential calcula- 
tions have been determined every 20 (16 at ft = 0.1565) 
trajectories. The numbers of thermalised configurations 
analysed in total are denoted by n g i uc and n pot in the 
Table for the two classes of measurements, respectively. 

In order to reduce statistical noise we had to employ 
measurement frequencies that are larger than the inverse 
integrated autocorrelation times in most cases. This is 
particularly true for the glueballs and torelons. Autocor- 
relation effects have been taken care of by binning the 
time series into blocks prior to the statistical analysis. 
The bin sizes have been increased until the statistical er- 
rors of the fitted parameters stabilised. The numbers of 
effectively independent configurations that, after this av- 
eraging process, finally entered the analysis are denoted 
by m g iue and m po t, respectively. 

In addition to the dynamical quark simulations, 
quenched reference potential measurements have been 
performed. For this purpose, smeared Wilson loop data 
generated at (3 = 6.0 and (3 = 6.2 in the context of the 
study of Ref. p4j| have been re-analysed. The hybrid IL U 
potential reference data have been obtained in Ref. |35| ] 
at the same two (3 values. 

TABLE I. Simulation parameters. The last two rows refer 
to quenched simulations. 



K 


LlL T 


r a 1 


m^aLcr 






flpot 


nipot 


0.1560 


16 3 32 


5.11(3) 


7.14(4) 


2128 


266 


236 


236 


0.1565 


16 3 32 


5.28(5) 


6.39(6) 


500 


250 


322 


161 


0.1570 


16 3 32 


5.48(7) 


5.51(4) 


1000 


200 


250 


125 


0.1575 


16 3 32 


5.96(8) 


4.50(5) 


2272 


142 


270 


90 


0.1575 


24 3 40 


5.89(3) 


6.65(6) 


1980 


110 


201 


67 


0.1580 


24 3 40 


6.23(6) 


4.77(7) 


1780 


89 


196 


49 


/3 = 6.0 


16 3 32 


5.33(3) 








570 


570 


13 = 6.2 


32 4 


7.29(4) 








116 


116 
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In order to determine glueball masses, temporal cor- 
relation functions between linear combinations of Wilson 
loops have been constructed. For simplicity we restricted 
ourselves to plaquette-like operators that can be built 
from four (fat) links. Torclons, i.e. flux loops that live on 
the torus and encircle a periodic spatial boundary [j36| , 
have been created from linear combinations of smeared 
spatial Polyakov (or Wilson) lines. Both, the fat pla- 
quettes, used to construct the glueballs, and the Wilson 
lines can have three orthogonal spatial orientations. The 
real part of the sum over these orientations transforms 
according to the + representation |37j of the relevant 
cubic point group, Oh <8> Z%. In order to project out the 
T± representation one can take the imaginary part of 
one of the three possible orientations; the E ++ represen- 
tation can be achieved by taking the real part, either of 
the sum of two orientations minus twice the third one or 
of the difference between two orientations. In order to 
increase statistics we average over the correlation func- 
tions obtained from the three equivalent orthogonal 
or two E ++ realisations. Within all three channels, + , 
E ++ and Tj , we subtract the disconnected parts from 
the correlation functions. This is only necessary for the 
A~l + state that carries the vacuum quantum numbers. 
However, in doing this for all channels, we find (slightly) 
reduced statistical fluctuations. 

All operators are projected onto zero momentum at 
source and sink by averaging over all spatial sites. In 
the case of glueballs, the Af + , E++ and T}~ represen- 
tations can be subduced from the continuum 0(3) ® 
J PC = ++ ,2 ++ and l" 1 representations, respectively, 
while the cubic point group remains the relevant symme- 
try group for the torelons (that only exist on a torus), 
even in the continuum limit. In what follows we shall 
label glueballs by the above continuum J PC quantum 
numbers. 

The fat links that form the basis of the glueball and 
torelon creation operators have been constructed by al- 
ternating "APE smearing" |3^] (keeping the length of 
the smeared link constant) with "Teper fuzzing" |3{| (in- 
creasing the length by a factor two). In both algorithms 
a link of smearing level i + 1 is formed by taking a linear 
combination of its predecessor of iteration i (in the case 
of fuzzing the product of two straight i level links) and 
the four closest surrounding spatial staples, followed by 
a projection back into the gauge group |Q. We denote 
the un-smeared links by i — 0. The first iteration then 
consists of APE smearing with the relative weight of the 
straight line connection with respect to the staple sum 
being a — 0. This is followed by a fuzzing iteration with 
a = 1.3, APE smearing with a = and so forth. This 
procedure is the outcome of a semi-systematic study of 
different fuzzing algorithms . Fat links of levels and 
1 have effective length a, levels 2 and 3 length 2a, levels 
4 and 5 length 4a and levels 6 and 7 length 8a. Itera- 
tions 8 and 9, yielding fat links of extent 16a, are only 
performed on the 24 3 lattices. 

APE smearing has also been applied to obtain ground 



state and excited state potentials. In this case, we iterate 
the procedure 26 times with a straight line coefficient, 
a = 2.3. Subsequently, Wilson loops are constructed 
whose spatial parts are built from smeared links. In ad- 
dition to on-axis distances, r || (1,0,0), off-axis sepa- 
rations that are multiples of (1,1,0), (1,1,1), (2,1,0), 
(2,1,1) and (2,2,1) have been realised for the ground 
state potential. In the case of the hybrid potential, in 
addition to on-axis separations, the plane-diagonal di- 
rection, r |j (1,1,0), has been investigated. The spatial 
gauge transporter along the latter direction is obtained 
by calculating the difference between the positively and 
negatively oriented "one-corner" paths connecting the 
two end points. For on-axis separations the potential in 
the continuum H u representation can be obtained from 
the E u representation of D^h- The E u creation operator 
has been constructed in the usual way |40[| from the sum 
of two differences between forwards and backwards ori- 
ented staples. The orthogonal depth of these staples has 
been chosen to be one lattice unit for r — a, two lattice 
units for 2a < r < 5a and three lattice units for r > 6a. 



III. MASS DETERMINATIONS 

A cross correlation matrix between the eight (ten in 
the case of the L a = 24 glueballs) basis operators cor- 
responding to different smearing/fuzzing levels is formed 
for each temporal separation, t. For the sake of numerical 
stability we restrict our analysis to two by two subsets^ of 
the correlation matrix, M(i), with blocking level i > 2. 
We then select the subset of two states for which the 
combination, M _1 / 2 (a)M(2a)M~ 1 / 2 (a), yields the low- 
est energy eigenvalue, and trace the t dependence of the 
effective mass, 



meff(t) 



111 



C(t) 
C(t + c 



(2) 



calculated from the correlation function, C(t), that is as- 
sociated with the optimal eigenvector, in order to identify 
a plateau. The glueball and torelon masses are subse- 
quently determined by means of a fully correlated and 
bootstrapped fit to the plateau region. The required 
covariance matrices are determined by means of sub- 
bootstraps. 

We denote the state that is created by an operator 
of fuzzing level i — 0, 1 , . . . within each glueball channel 
by \i). In Tabs. || and III the state vectors used to ex- 
tract the ++ and 2 ++ glueballs are listed, respectively, 
as well as the fitted overlaps of these states with the phys- 
ical ground states. We have employed the normalisation, 
= 1. The data at K = 0.1565 is less precise due to the 



In the case of the 1 + glueball only the diagonal elements 
of the correlation matrix are used. 
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smaller ensemble size, n g i uc . In general, our variational 
procedure yields ground state overlaps around 80 %. 

While a direct determination of glueball radii turns 
out to be impossible within our limited statistical en- 
semble sizes, the shapes of the optimal creation opera- 
tors are indicative for the approximate diameters of the 
underlying physical ground states. Note that the only 
direct measurement of glueball sizes to-date has been at- 
tempted in pure SU(2) gauge theoryin Ref. ^d|. On all 
16 3 volumes we find the ++ wave function to receive 
a dominant |5) component, which is one APE iteration 
applied to a fuzzed link of effective length 4a. The sub- 
leading contribution is |6), i.e. a square with side length 
8a. We conclude that the "diameter" of the scalar glue- 
ball is somewhat less than 8 lattice spacings or 0.7 fm. 
These observations indicate a somewhat wider glueball 
than observed in previous quenched studies Hl4jj,^l],0] . 

As one increases k, the effective lattice spacing, deter- 
mined for instance in units of the gluonic observable, r , 
decreases and the glueball extends over more lattice sites. 
This results in an increased relative weight of the larger, 
1 6), state, which finally becomes the dominant contribu- 
tion at k = 0.158. There is a difference between the 
two volumes at k = 0.1575: the size of the glueball on 
the larger lattice becomes smaller, indicating finite size 
effects. As we shall see below this goes along with an 
increase of its mass. 

The 2 ++ state is dominated by the |6) contribution, 
with the exception of k = 0.1565. At this k value the 
statistics are not really sufficient for a reliable determina- 
tion of the optimal state vector, as indicated by the large 
relative error on the ground state overlap, 0.48 ± 0.20. 
The second largest contribution then comes from the |4) 
(|3) for the k = 0.1575, L a = 24 simulation) state. In 
general, we find the |4) state to have best overlap with 
the l" 1 glueball. 



TABLE II. 


Scalar 


glueball wave function and 


overlap with 


the physical ground state. 






L a 


state vector 


g.s. overlap 


0.1560 


16 


0.76|5) +0.65|6) 


0.89(2) 


0.1565 


16 


0.77 5} +0.63 6} 


0.71(12) 


0.1570 


16 


0.72 5} +0.70 6} 


0.79(4) 


0.1575 


16 


0.73 5} +0.68 6} 


0.83(5) 


0.1575 


24 


0.72 4} +0.70 5} 


0.94(9) 


0.1580 


24 


0.65 5} +0.76 6} 


0.80(8) 
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FIG. 1. Effective scalar and tensor glueball masses and fits 
at k = 0.156. 
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FIG. 2. Same as Fig. for k = 0.1565. 



TABLE III. Tensor glueball wave function and overlap 
with the physical ground state. 



K 


L a 


state vector 


g.s. overlap 


0.1560 


16 


0.37|4} +0.93|6) 


0.81(3) 


0.1565 


16 


0.83 4} +0.55 6} 


0.48(20) 


0.1570 


16 


0.28 4} +0.96 6) 


0.81(3) 


0.1575 


16 


0.27 4} +0.96 6} 


0.80(3) 


0.1575 


24 


0.37 3} +0.93 6} 


0.62(17) 


0.1580 


24 


0.66 4} +0.75 6} 


0.84(4) 



4 
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FIG. 3. Same as Fig. g for k = 0.157. 
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Fig 


. for the 16 3 


volume at k — 


0.1575. 




TABLE IV. Torelon ground state overlaps 




K 


La 






E++ 


r+- 


0.1560 


16 




0.89(2) 


0.90(2) 


0.91(2) 


0.1570 


16 




0.87(6) 


0.82(4) 


0.72(11) 


0.1575 


16 




0.41(10) 


0.74(2) 


0.39(12) 


0.1575 


24 




0.55(13) 


0.84(3) 


0.60(21) 


0.1580 


24 




0.85(24) 


0.83(11) 


0.24(9) 



In Figs. [I] - |[ we show effective mass plots for the 
scalar and tensor glueballs at the various parameter val- 
ues, together with the corresponding fit results (solid 
lines with dashed error bands). Plateaus could be identi- 
fied either from t = 2a or t = 3a onwards. Since we take 
correlations between the data points into account fitted 
results can turn out to lie somewhat below the central val- 
ues of the individual effective masses (cf. Figs. | and |). 
The torelon masses have been analysed in an analogous 
way. All torelons receive their dominant contribution ei- 
ther from the |6) or |7) channels with a sublcading term 
from the |4) or |5) states. The achieved ground state 
overlaps are listed in Tab. |y|. The statistical quality of 
the signal did not permit an estimation of the overlap at 
k = 0.1565. 
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FIG. 5. Same as Fig. |l| for the 24 3 volume at K = 0.1575. 
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FIG. 6. Same as Fig. [l] for k = 0.158. 

On volumes of similar physical size in terms of r$ we 
find statistical errors of fuzzed plaquettes and Wilson 
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lines to be smaller than in comparable quenched simu- 
lations P^lfl- This can very well be due to the larger 
number of degrees of freedom when including quarks. 
Similar conclusions can be drawn from a comparison be- 
tween published SU(2) and SU(3) results. In the case of 
the determination of the static potential absolute errors 
on smeared Wilson loops are smaller than in comparable 
quenched simulations too. However, this effect is partly 
compensated for by a faster decay of the signal in Eu- 
clidean time, due to an increase of the static source self 
energy. In most gluonic quantities self-averaging over 
more lattice points when increasing the physical volume 
reduces the statistical noise. This is of course not the case 
for torelons that become heavier with growing volume. 
However, statistical fluctuations of glueball correlation 
functions grow with the volume too. This effect, that 
has also been observed in quenched simulations 0,0, 
is related to contributions stemming from large distance 
correlators to the spatial sum that is required to project 
onto zero momentum. 

In the potential measurements no variational optimi- 
sation is performed. Again, effective masses are traced 
against t, separately for all lattice separations, r. In view 
of the fact that, given the exponential increase of rela- 
tive errors with t, the result of a fit to data at t > t m i n 
will be almost identical to the effective mass calculated 
at t — t m in, we decided to estimate the potential by the 
latter. We also demanded the effective mass at t m i n — a 
to be compatible with the quoted result, determined at 

t tmin ■ 



In Fig. 
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we illustrate the convergence of effective 
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V(r,t) = a -1 In 



W(r,t) 
W(r,t + a) 



V(r) (t-oo), (3) 



towards the asymptotic value, V(r), for the case of the 
re = 0.1575 potential on the 24 3 volume at four randomly 
selected distances. Note that, due to the positivity of the 
Wilson action, the approach towards the plateau has to 
be monotonous. At (3 — 5.6, with the smearing algorithm 
used, the range of re values investigated and the selection 
criterium employed, we found V\r, t m i n (r)] to approxi- 
mate the asymptotic value within statistical errors for 
4a < tmin < 6a, i.e. at i m m ~ fg. In the case of the TL U 
potentials, the ground state overlaps were found to be 
somewhat inferior and the statistical errors larger. As a 
consequence of these two competing effects the necessary 
t-cuts came out to be almost identical to those employed 
for the ground state potential. 

For distances bigger than about 2ro, t m in bad to be 
increased by one lattice unit with respect to its value at 
smaller distances. This differs from our experience from 
quenched studies employing similar methods at similar 
lattice spacings where systematic and statistical effects 
happened to interfere in such a way that the same t-cut 
could be applied at all distances 13 lj] . Once the po- 
tential is extracted, the overlap of the creation operator 
with the physical ground state can be quantified in terms 
of an overlap coefficient, < Co(r) < 1: 

C (r,t) = W(r,t)e v ^ t ^ C (r) [t -» oo). (4) 



o 




FIG. 7. Effective masses for the potential on the 24 3 lattice 
at k= 0.1575 at four different source separations. 



FIG. 8. Ground state overlaps of the potential creation op- 
erators as a function of r. The quenched results have been 
obtained at (3 = 6.2. 

In Fig. ||, we compare the ground state overlaps that 
have been achieved at re = 0.156 and re = 0.1575 
(La = 24) with quenched (3 = 6.2 data that have been 
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obtained by use of the same smearing algorithm and anal- 
ysis procedure. We note that while at small r the ground 
state overlaps of the un-quenched simulations are supe- 
rior to the quenched reference data, at large r the projec- 
tion onto the physical ground state becomes increasingly 
worse, in particular at the lighter quark mass. 

The dependence of the quenched overlaps on r is lin- 
ear to first approximation, however, this is not so for 
the sea quark data (which explains why t m - ln had to be 
increased as a function of r in the latter case) . We con- 
clude that at large r the flux tube distribution appears 
to change when including sea quarks. This effect, which 
becomes more pronounced when the quark mass is de- 
creased, can be interpreted as a first indication of string 
breaking: the physical ground state will be a mixture 
between a would-be static-static state and a pair of two 
static-light mesons. While our creation operator is tuned 
to optimally project onto the former it has almost zero 
overlap with a wave function of the latter type. As a 
consequence, we lose local mass plateaus altogether in 
the regime, r > 2.5 ro. A similar reduction of ground 
state overlaps at large distances has been reported by 
the CP-PACS Collaboration @ . 

IV. RESULTS 



r > 2ro, along with a study of the corresponding tran- 
sition matrix elements that might shed more light onto 
this question @|^. 
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0.5 1 1.5 2 2.5 3 

r/r 

FIG. 9. Comparison between quenched (/3 = 6.2) and 
un-quenched ground state and il„ potentials. 



A. The static potentials 

The potentials obtained from the six dynamical 
fcrmion and two quenched simulations have been rescaled 
in units of ro. Subsequently, the constant roV(ro) has 
been subtracted in all cases to cancel the (different) static 
self energy contributions and to achieve the common nor- 
malisation, V(tq) = 0. Up to violations of rotational 
symmetry at small distances, the un-quenched potentials 
are found to agree with each other within statistical er- 
rors. In particular, we do not see finite size effects among 
the two k = 0.1575 data sets. 

In Fig. we compare the potential obtained on the 
largest lattice volume at our disposal, L a a = 24 a « 
4.07 ro « 2.03 fin at k = 0.1575, with the quenched po- 
tential at (3 — 6.2. Besides the ground state potentials, 
which correspond to the representation of the rele- 
vant continuum symmetry group -Doo/i, the II U hybrid 
potentials jl0],U are depicted as well as the approximate 
masses of pairs of static-light scalar and pseudoscalar 
bound states into which the static-static string states 
are expected to decay at large r. Note that for sep- 
arations, r > 1.9 ro, we have used some extra off-axis 
directions, in addition to those mentioned in Sect. 0. 
Around r w 2.3 ro we expect both, and Tl u potentials 
to exhibit colour screening. However, the Wilson loop 
data are not yet precise enough to resolve this effect. In 
a forthcoming paper p5[ we therefore intend to present 
a more systematic analysis geared to improve the statis- 
tical quality of potential data in the interesting regime, 
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r/r 

FIG. 10. First hybrid excitation: quenched vs. un-quenched. 

In Fig. |l^, we compare our two most precise un- 
quenched data sets (k = 0.156 and k = 0.1575) on the 
hybrid potential with quenched results. We find no sta- 
tistically significant differences, other than finite a effects 
between the quenched potentials obtained at two differ- 
ent (3 values (full symbols). 

The ground state potentials are fitted to the parametri- 
sation, 



V(r) = V + ar- 
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where the term [1/r], that denotes the tree level lattice 
propagator in position space [[|, is included to quantify 
the short distance lattice artefacts. Given the fact that 
the number of elements in the covariance matrix is much 
bigger than the sizes of our data samples, we refrain from 
attempting correlated fits. Therefore, x 2 can in principle 
become much smaller than the degrees of freedom, Ndf- 
All errors are bootstrapped. 

We wish to investigate the effect of sea quarks onto the 
short range interquark potential. Apart from the fit pa- 
rameters, we also calculate the Sommer scale of Eq. (pi), 



1.65 



(6) 



The fit function, Eq. (||) , is only thought to effectively pa- 
rameterise the potential within a given distance regime, 
rather than being theoretically sound. For instance nei- 
ther string breaking nor the running of the QCD coupling 
has been incorporated. Therefore, the values of the effec- 
tive fit parameters will in general depend on the fit range, 
[fiain, r roax ]. For instance, as a consequence of asymp- 
totic freedom, e will weaken as data points from smaller 
and smaller distances are included. To exclude such a 
systematic bias from our comparison between quenched 
and un-quenched fit parameters, we deviate somewhat 
from our first analysis (as we have already done in 
Refs. [pp|^l| ) and use the same r min > 0.4 r in physical 
units for all our fits and r max < aL a /2. The results are 
displayed in Tab. [v| 

Three-parameter fits have been performed in addition, 
by constraining g = in Eq. (||) . To obtain acceptable fit 
qualities for the same physical fit range on all data sets 
we had to increase, r m i n ps 0.6 tq. The latter fit results 



are shown in Tab. VI 



TABLE V. Four-parameter 
r m in > 0.4 ro- 



fits to the static potential for 



0.1560 
0.1565 
0.1570 
0.1575 
0.1575 
0.1580 



La 


r min 
a 


x 2 


9 


V a 


16 


V5 


0.94 


0.27(5) 


0.731(06) 


16 


^3 


0.54 


0.42(7) 


0.738(10) 


16 


75 


0.78 


0.34(7) 


0.725(12) 


16 


Ve> 


0.47 


0.30(9) 


0.739(11) 


24 




0.76 


0.30(6) 


0.746(05) 


24 


2^2 


0.64 


0.28(7) 


0.751(07) 










0.760(12) 



0.338(09) 
0.349(16) 
0.324(19) 
0.337(20) 
0.352(09) 
0.358(12) 
0.368 



(3 = 6.0 16 
(3 = 6.2 32 



3 



+ 20 
26 



0.224(2) 
0.216(3) 
0.210(4) 
0.192(4) 
0.193(2) 
0.182(3) 
0.171 



+6 
3 



0.14 0.34(2) 
0.49 0.25(1) 



0.658(04) 
0.636(03) 



0.291(05) 
0.292(05) 



0.219(2) 
0.160(1) 



In Tab. VIl, we display the ro values obtained from 
our four-parameter interpolations as well as the effective 
string tension in units of ro- The same data is depicted 
in Fig. ^ as a function of the squared it mass, (m % a) 2 . 
Chiral extrapolations are performed according to, 

ar Q _1 (m 7r ) = ar^ 1 (0) + cia 2 + c 2 a 4 m^. (7) 

The fit to all data points, with the exception of the L a — 
16 result at k = 0.1575, yields the parameter values, 



aro X (0) = 0.1476(30), 
ci = 0.324(57), 
c 2 = -0.41(22), 

X 2 /N DF = 1.37/2. 



(8) 
(9) 
(10) 

(11) 



When excluding the rightmost point [k = 0.156) a linear 
extrapolation becomes possible: 



ar-^O) = 0.1509(14), 

ci = 0.246(16), 
X 2 /N DF = 1.32/2. 



(12) 
(13) 
(14) 



The fitted curves as well as the quadratically extrapo- 
lated Tq 1 (0) value (open square), with error bars enlarged 
to incorporate the one a region around the linear extrap- 
olation, are included in the Figure. 



TABLE VI. 
nm « 0.6 ro. 



Three-parameter fits to the static potential for 



K 




r min 
a 


x 2 


Voa 


e 




0.1560 


16 


3 


0.84 


0.712(12) 


0.290(23) 


0.228(3) 


0.1565 


16 


3 


0.48 


0.715(21) 


0.291(42) 


0.221(6) 


0.1570 


16 


3 


1.43 


0.725(21) 


0.326(42) 


0.210(6) 


0.1575 


16 


3 


0.57 


0.727(19) 


0.308(38) 


0.195(5) 


0.1575 


24 


3 


0.75 


0.736(07) 


0.322(15) 


0.195(2) 


0.1580 


24 


3 


0.77 


0.747(08) 


0.348(17) 


0.183(3) 


/3 = 6.0 


16 


3 


0.87 


0.639(06) 


0.245(19) 


0.223(3) 


(3 = 6.2 


32 


3V2 


0.67 


0.637(05) 


0.296(13) 


0.160(1) 



TABLE VII. r and the ratio ^/ar . 



K 


La 


ro/a 


s/ffVQ 


0.1560 
0.1565 
0.1570 


16 
16 
16 


5.104(29) 
5.283(52) 
5.475(72) 


1.145(4) 
1.141(7) 
1.152(8) 


0.1575 
0.1575 
0.1580 


16 

24 
24 


5.959(77) 
5.892(27) 
6.230(60) 
6.73llg 


1.146(9) 
1.139(4) 
1.137(5) 
1.1331" 


13 = 6.0 
13 = 6.2 


16 
32 


5.328(31) 
7.290(34) 


1.166(2) 
1.165(2) 
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FIG. 11. Chiral extrapolation of the scale r^" 1 . The open 
square corresponds to the chiral limit, the full square to the 
physical ratio, rom^ w 0.350. 



An extrapolation to the physical limit, m n h ro 
0.350, yields, 



ar - 1 (m 7rp J = 0.1485^. 



(15) 



The corresponding r$/a value is included in Tab. VII. 
Note that n p h — 0.158456(19) only marginally differs 
from k c = 0.158493(18). 

From bottomonium phenomenology [ppl|j34[| one ob- 



tains, r. 







(394 ± 20) MeV. The lattice spacing de- 
termined from tq at physical sea quark mass, therefore, 
is, 



(2.68±? x ± 0.14) GeV. 



(16) 



The last error reflects the scale uncertainty within the 
phenomenological Tq determination. The above value 
compares well with a -1 = 2.7(2) GeV, as obtained from 
the p mass after an extrapolation to the phys- 

ical m^/vrip ratio. By interpolating between quenched 
reference data ||, we find that the above tq value cor- 
responds to the quenched [3 = 6.14(2): inclusion of two 
light Wilson quarks at j3 = 5.6 thus results in a j3 shift, 
A/3 = 0.54(2), i.e. in an increase of the lattice coupling, 
g 2 = 6/(8, by 9-10 %. 

The values of Vga, y/aa and e, extrapolated to the 
physical quark mass, are displayed in Tab . fv| while the 
combination y/aro is included in Tab. [VI]] . In Figs. [l2| 
and [b| we display y/aro and e against k, respectively. 
In these two cases not only quadratic and linear fits in 
(m w a) 2 but also simple averages have been taken into ac- 
count when assigning the error bars, that reflect both sta- 
tistical errors as well as the uncertainty in the parametri- 
sation, to the (quadratically) extrapolated values. 



1.17 
1.16 
1.15 
1.14 
1.13 
1.12 



quenched 0.156 



0.157 



0.158 



FIG. 12. Sea quark mass dependence of the ratio between 
the square root of the string tension and r^ 1 . The horizontal 
lines with error bands correspond to the quenched (solid) and 
chirally extrapolated un-quenched (dashed) results. 
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FIG. 13. Sea quark mass dependence of the effective 
Coulomb strength of the static potential. 



From Fig. [L2] as well as from Tab. VII we can read 



off that the ratio y/aro decreases down to 1.133I 8 , com- 
pared to the quenched result, 1.165(3). With r w 0.5 fm 
we obtain a quenched value, yfa 460 MeV, while with 
two sea quarks we find, y/a « 445 MeV. Note that the lat- 
ter result comes closer to estimates, Ja — (429±2) MeV, 
from the p, ct2, . . . Regge trajectory ||. 

We find self energy, Vba, and effective Coulomb co- 
efficient, e, to significantly increase with respect to the 
quenched simulations. In tree level perturbation theory 
both parameters, V^a and e, are expected to be propor- 
tional to g 2 . Data on e are shown in Fig. [l3| The re- 
sults, extrapolated to the physical point, are (Tab. 0), 
e = 0. 368^26 an d Vqo, — 0.760(20), i.e. e is increased 
by 16 - 33 % and Vq by 16 - 21 % with respect to the 
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quenched results: the effect of sea quarks on couplings de- 
fined through the potential at short range is bigger than 
the relative shift in the lattice coupling, A/3//3 0.1. 
Bottomonia spectroscopy suggests a value [Q, e ps 0.4, 
which, given our nf — 2 result, is indeed likely to be 
consistent with QCD with three light sea quark flavours. 
The three-parameter fits (Tab. |Vlj ) are, due to the differ- 
ent fit range r > 0.6 ro, not yet precise enough to confirm 
an increase of e beyond doubt. However, at least the ef- 
fect on Voa is statistically significant and in agreement 
with the result from the four-parameter fits. 

In Fig. |lj, we compare quenched (open symbols) and 
un-quenched (full symbols) lattice data at short dis- 
tances to visualise that our interpretation of an increase 
in the effective Coulomb strength is independent of the 
fitted parametrisation. The curves correspond to the 
values of the parameters org, e and Vorg at ft = 6.2 
(quenched) and k — 0.1575 (un-quenched) that have 
been obtained in the four-parameter fits. Due to lat- 
tice artefacts the data sets scatter around the interpo- 
lating curves. As expected, we indeed observe that the 
un-quenched data points (full symbols) lie systematically 
below their quenched counterparts (open symbols). 
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FIG. 14. The short distance potential: quenched (open 
symbols) vs. un-quenched (full symbols). 



B. Torelons 



in the limit of large volumes, aL a — > oo. a should be the 
very same string tension that governs the static poten- 
tial at large distances. The sub-leading finite size correc- 
tion il, 



Act = - 



3(aL CT ) 2 



(18) 



is expected from the bosonic string picture. The effect of 
this contribution, which has been accurately verified in 
SU(2) gauge theory p2[ , ranges from a decrease by 11 % 
on our smallest physical volume (16 3 at k = 0.1575) to 
4.9 % on our biggest volume (24 3 at k = 0.1575). 

In a pure gauge theory the centre symmetry of the 
action implies that all torelons that wind once around a 
lattice boundary are mass degenerate (see e.g. Appendix 
D.2 of Ref. || for details). When explicitly violating the 
centre symmetry by including sea quarks the degeneracy 
will be broken, with torelons to be classified in accord 
to representations of the cubic group, Oh times charge 
conjugation. While in pure gauge theories glueballs can 
only split up into pairs of torelons, with sea quarks a 
single torelon can mix with or decay into a glueball that 
is in the same representation. For this to happen the 
torelon has to break up, unwind and rearrange itself into 
a flux loop with trivial winding number. In reality, one 
will meet different states (and creation operators) with 
equal quantum numbers. The operator that corresponds 
to the torelon in the quenched case will have maximal 
overlap with the state that is most torelon like but will 
in general, in the limit of large Euclidean times, decay 
into the lightest possible state. In the case that another 
state turns out to be lighter than the expected torelon 
mass we would call this effect "string breaking" of the 
torelon. 

Heuristically one would expect a torelon-type opera- 
tor to project better onto a state with a large mesonic 
component than a state with a large glueball component 
since the former appears as an intermediate state in the 
deformation of a torelon into a glueball anyway. On the 
other hand progress is being made in the excitation of 
flavour singlet states by use of mesonic operators (quark 
loops) (5^|l(],|54| and there are indications |l(J that purely 
gluonic operators of the type we use as well as mesonic 
operators might both have acceptable overlaps with the 
physical ground state. 



In view of the difficulties in detecting string breaking 
directly in the static potential, it might be illuminating 
to investigate the effect of including sea quarks on torelon 
masses pqj . Such numerical simulations have been pio- 
neered by Kripfganz and Michael on small lattice 
volumes some ten years ago. In pure gauge theories a 
torelon that winds around a dimension of extent L a a 
will have a mass, 



TABLE VIII. Torelon masses. 



L a m.++r m B+ +r m T+ -r expected 2m n ro 



0.1560 
0.1565 
0.1570 
0.1575 
0.1575 
0.1580 



16 
16 
16 
16 
24 
24 



3.98(30) 
3.34(36) 
3.53(25) 
2.85(17) 
4.41(68) 
6.11(93) 



4.03(20) 
3.43(25) 
3.54(19) 
2.82(15) 
6.16(66) 
5.36(44) 



3.84(32) 
3.59(42) 
3.06(33) 
2.82(18) 
4.06(63) 
2.24(98) 



3.78(07) 
3.59(11) 
3.51(14) 
3.11(12) 
5.03(08) 
4.71(13) 



4.56(4) 
4.22(6) 
3.77(6) 
3.32(5) 
3.26(3) 
2.47(5) 



nrix — aL a [a + Aa(aL a 



(17) 
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FIG. 15. Effective torelon masses at k = 0.156, in compar- 
ison with the expectation, Eqs. (|l7]) - (|l8|), with the effective 
string tension cr obtained from a fit to the static potential. 
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FIG. 16. Same as Fig. |l|at « = 0.1575. 

In addition to the fitted torelon masses, we display the 
non-s tring breaking expectations of Eqs. ( |l7| ) - ( [l8| ) in 
Tab. VIII, based on the string tensions from our four- 
parameter fits to the static potential of Tab. Unfortu- 
nately, the signals on the large lattices become very noisy, 
due to the larger torelon mass. Apart from isoscalar me- 
son/glueball states, the Af + torelon can in principle de- 
cay into a pair of 7r's. We therefore include for guid- 
ance twice the tt mass in the last column of the Table. 
As we shall see below (Tab. IX), on all 16 3 lattices, the 



scalar glueball comes out to have a mass very close to 
the torelon mass expectation of Eqs. (17) and (|l8|). On 
our 24 3 lattices, however, the expectation supercedes the 
++ and, eventually, the l" 1 glueball masses: breaking 
of the and, eventually, torelons into a flavour 
singlet meson/glueball state becomes energetically pos- 
sible. At the same time two 7r's drop below the torelon 



k = 0.1575, L„ = 24 



expectation (and become lighter than the flavour singlet 
meson/glueball state). 
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FIG. 17. Same as Fig. |l^ on the larger volume. 
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k = 0.1580, L„ = 24 
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FIG. 18. Same as Fig. at k = 0.158. 

On all the 16 3 lattices, within statistical errors, the 
three torelons are mass degenerate and in agreement with 
the non-string breaking expectation. Two effective mass 
plots are displayed as examples in Figs. [l5| and |l^. The 
same holds true, within larger errors, in the case of the 
24 3 simulation at k = 0.1575 (Fig. |l7|). However, at 
k = 0.158 with a mass of (2.24 ± 0.98) r at least the 
torelon appears to fall short of the expectation by 
two and a half standard deviations (Fig. [L8|)! It would 
be important to increase statistics in order to corroborate 
this result as a signal for string breaking. 
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C. Glueballs 

We display the extracted masses of the three glueballs 
that we investigated in Tab. |x|. The quenched results 
at a similar lattice spacing (/3 = 6.0) have been calcu- 
lated by Michael and TeperQ [0 . The quenched contin- 
uum 0++ result § (last row of the Table) stems from 
a quadratic extrapolation in the lattice spacing of data 
obtained in Refs. |55| , ^6| , ^2| j2U5^] while the tensor and ax- 
ialvector results are from Morningstar and Peardon [Q. 
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FIG. 19. Effective 1 glueball masses. The horizontal line 
corresponds to the quenched continuum result. 
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FIG. 20. Results on scalar and tensor glueball masses. The 
horizontal lines with error bands indicate quenched contin- 
uum results while the leftmost data points correspond to 
quenched Wilson action results obtained at finite lattice spac- 
ing (f3 = 6.0). The expected torelon masses on 16 3 and 24 3 
lattices are shown too (open triangles). 

Effective mass plots of the scalar and tensor channels 
are shown in Figs, [l] - ^. In the case of J PC = R no 
plateau has been established and the values displayed in 
the Table only represent upper limits. The corresponding 
effective masses are displayed in Fig. [r], in comparison to 
the continuum limit extrapolated quenched result (hori- 
zontal line with error band). Quenched results obtained 
at similar lattice spacings jpj suggest the glueball to lie 
somewhere inbetween 6.5 and 8.5 Tq 1 and, within sta- 
tistical errors, do not deviate from the continuum value. 
Keeping in mind that all un-quenched simulations are 
statistically independent of each other, there is a clear 
indication for a R state being lighter than 6r ~ 1 . It re- 
mains to be demonstrated in future simulations at differ- 
ent lattice couplings, (3, whether this conclusion persists 
in the continuum limit. 

In Fig. our results on the scalar and tensor glueballs 
are summarised. The leftmost data points represent the 
(3 = 6.0 quenched results while the horizontal error bands 
correspond to the respective quenched continuum limits. 
The quenched scalar glueball at the finite lattice spacing 
comes out to be lighter by 15-20 % than the continuum 
limit extrapolated value |^,|5| while the tensor glueball 
is in agreement with the continuum result. In addition 
to the glueball masses, the torelon mass expectations of 
Eqs. ( p7| ) - fll8| ) are included into the Figure. For k < 
0.157, 2m„ is larger than the torelon mass. On the 16 3 
lattice at n = 0.1575 the two masses come out to be 
degenerate within errors and on both 24 3 lattices two 7r's 
become lighter than the torelon (and the scalar glueball). 

TABLE IX. Glueball masses. A star indicates that no ef- 
fective mass plateau could be identified and, therefore, only 
an upper limit is stated. 



ft 




m ++r 


m 2++ r 


m 1+ _r 


m 2++ /m 0++ 


0.1560 


16 


3.56(12) 


5.77(19) 


5.0(2.1)* 


1.62(07) 


0.1565 


16 


3.72(46) 


4.97(106) 


8.1(1.2)* 


1.34(34) 


0.1570 


16 


3.02(18) 


5.70(19) 


4.6(1.2)* 


1.89(13) 


0.1575 


16 


3.26(25) 


6.06(23) 


12.7(2.3)* 


1.86(13) 


0.1575 


24 


4.21(35) 


5.90(82) 


6.2(2.3)* 


1.41(23) 


0.1580 


24 


3.85(38) 


6.65(36) 


4.8(1.8)* 


1.72(17) 


,5 = 6.0 


16 


3.57(21) 


5.70(32) 


7.5(1.1) 


1.60(14) 


/3 = 00 




4.36(7) 


5.85(8) 


7.18(8) 


1.34(3) 



2 In the case of the 2 ++ glueball we display their E ++ result 
only since the T 2 + representation has not been realised in 
our study. 
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We do not detect any statistically significant deviation 
between quenched and un-quenched results in the 2 ++ 
channel. Unfortunately, the statistical error on the 24 3 
lattice at k — 0.1575 is too large to resolve possible finite 
size effects. We find all 16 3 ++ glueball data to roughly 
agree with the quenched (3 — 6.0 result. At the same 
time, within statistical errors, the scalar glueball masses 
are degenerate with the Af + (expected and measured) 
torelon masses. This degeneracy has also been observed 
in Ref. |50| , on small lattice volumes: in the quenched 
approximation, finite size effects are small as long as two 
torelons are heavier than the glueball, however, with sea 
quarks, there is no protection preventing the glueball to 
decay into a single torelon. This could explain why the 
++ glueball on the 2 fin lattice at k = 0.1575 comes out 
to be significantly heavier than the one on the 1.4 fm lat- 
tice. Note that the size of the corresponding wave func- 
tion becomes reduced when increasing the lattice volume 
(Sect. 

On both 24 3 lattices the mass of the scalar glueball 
comes out to be bigger than 2m T . There is no physical 
reason why the creation operator employed should have 
zero overlap with either a pair of 7r's or a hypothetical tttt 
bound state |)8|. However, it seems that this overlap is 
small, similar to the situation of the static potential, in 
which the smeared Wilson loop only receives a tiny con- 
tribution from pairs of static- light states. In conclusion, 
the mass of the state onto which the scalar glueball op- 
erator dominantly projects is compatible with quenched 
glueball results obtained at j3 > 6. We find finite size 
effects that can be interpreted as mixing with torelon 
states. 



V. SUMMARY AND DISCUSSION 

We have determined various "pure gauge" quantities 
in a simulation of QCD. Including two sea quark flavours 
results in a shift of about 10 % in (3 and in an increase of 
the effective Coulomb strength, governing the interquark 
potential at r > 0.2 fm, by 16 - 33 %, in the limit of light 
quark masses. While we have not been able to establish 
string breaking either in the ground state or in the IT U 
hybrid potential, a comparison of the ground state over- 
laps with quenched reference data suggests an effect at 
large distances. 

Around r c w 2.3 ro w 1.15 fm string breaking will be- 
come energetically possible for both potentials, ground 
state and hybrid, at k, = 0.1575. On the assumption that 
reducing the bare quark mass, m = (k _1 — k~ )/(2a), by 
an amount, Am, induces an equal change in the static- 
light meson mass, we would guesstimate the string break- 
ing length scale to drop by Ar ss 2Am/a 0.18 ro when 
extrapolating to the physical sea quark mass. This crude 
estimate tells us that with two light quark flavours the 
QCD potential should become flat around r c ss 1.05 fm, 
such that in "real" 2 + 1 flavour QCD string breaking at 



distances of 1 fm or smaller is likely. 

For k < 0.1575 torelon masses are in agreement with 
the (non-string breaking) expectation in terms of the ef- 
fective string tension, obtained from a fit to the static 
potential, while at n — 0.158 there are indications for 
the torelon in one particular representation to become 
lighter than one would naively have expected. The basic 
problem is that we are only faced with a small window of 
lattice volumes for observation of "string breaking" . On 
a large lattice torelons become so heavy that the corre- 
lation functions disappear into noise at small temporal 
separation while on a small lattice the expected mass 
does not yet sufficiently differ from the masses of decay 
candidates to resolve an effect. It would be helpful in this 
respect to have for instance k = 0.1575 data at L a = 20. 
The tensor glueball comes out to be heavier than the re- 
spective torelon state, i.e. in this case the two creation 
operators have little overlap with each other. 

On our largest lattice (k = 0.1575) that corresponds to 
m^/nip — 0.704(5) with extent, L a a ~ 4.07 ro, we find, 
TOq++ = 4.21(35) Tq , which is in agreement with the 
quenched continuum limit result, m ++ — 4.36(7) Tq 1 . 
The 2 ++ glueball masses also seem to agree with the 
quenched results while there are indications for a light 
glueball. On small lattices we find the scalar glue- 
ball and torelon to be degenerate and the glueball mass 
is underestimated. To exclude such finite size effects, the 
spatial lattice extent should be taken to be bigger than 
about 3.5 ro ~ 1.75 fm. 

Previous results with two flavours of Kogut-Susskind 
fermions have been obtained by the HEMCGC Collab- 
oration H at the somewhat more chiral values [j59| , 
m^/rrip = 0.658(5) and m^/rrip = 0.434(4) on lattices 
with similar resolution as in the present study, ro ~ 5a. 
Their lattice extents are ]Tq], L a a = 3.08(12) ro and L a = 
2.50(5) ro, respectively, and, despite the rather small vol- 
ume in the latter case and the extremely light quarks in 
the former case, the results ||[l8), m ++ — 4.47(43) r^ 1 
and m ++ = 4.17(30) r^ 1 , are compatible with ours. 

Recently, the UKQCD Collaboration has reported 
first results for the scalar glueball, obtained with the 
Sheikholeshlami-Wohlert action on lattices with the lin- 
ear dimension, L a a ps 3.4 ro, at chiralities within the 
range explored in our study, as indicated by the ra- 
tios |To| , m^/nrip w 0.71 and m^jvrip ps 0.67. How- 
ever, their lattices are quite coarse: ro ~ 3.44 a and ro ~ 
3.65 a, respectively. The values they state are remarkably 
small, namely m ++ = 1.6(2) r^ 1 and m Q ++ — 1.8(4) Tq 1 , 
i.e. they come out to be by a factor two smaller than our 
16 3 lattice results, indicating larger finite lattice spacing 
effects than seen in quenched simulations f|[Iol- This 
issue certainly deserves further study. 

At sufficiently light sea quark masses we expect the 
scalar glueball to mix with either flavour singlet mesonic 
states and eventually tttt molecules or to decay into tt 
pairs. Such scenarios should be studied, in un-quenched 
as well as quenched simulations PJ5q] . In view of the mix- 
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ing and decay channels that open up once sea quarks are 
switched on, it is certainly also worthwhile to investigate 
glueballs with other quantum numbers than those that 
we have studied so far. In particular the pseudoscalar 
should be an interesting candidate, given the large mass 
of the rj meson. In view of fragmentation models as well 
as achieving an understanding of the decay rates of T (45) 
and T(55) states into pairs of BB mesons, a determi- 
nation of the string breaking distance and the associ- 
ated mixing matrix elements is certainly equally exciting. 
Work along this line is in progress fl6|,f45| . 
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